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Mr. W. D. Niven on the Calculation 


“ On the Calculation of the Trajectories of Shot.-” By W. 
D. Niven, M.A., F.R.A.S. Communicated by J. Clerk 
Maxwell, F.R.S., Professor of Experimental Physics in the 
University of Cambridge. Received March 24, 1876*. 

In the present state of our knowledge of the resistance of the air to 
shot, the problem of integrating the equations of motion of the shot and 
of plotting-out a representation of the curve described by it is peculiar, 
because, according to the best experiments we possess, the law of the 
retardation cannot be expressed by a single exact formula which is 
available for the solution. We are therefore compelled to give a 
solution adapted to Tables, the magnitudes of the retardation being set 
down in those Tables for velocities which are common in practice. The 
formulae given by Hutton and by Didion, even if they were true, apply 
only to spherical shot; and though they are very simple formula, the 
solutions obtained by means of them are not satisfactory—first, by 
reason of their complexity, and next on account of the rough approxima¬ 
tions which characterize the proofs. 

Prof. Helie, who gives an account of Didion’s method in his 4 Traite de 
Balistique/ says that it gives results which are not in accordance with 
fact. The fault may probably be laid in a great measure to the charge of 
the formula; for there can be no doubt that Mr. Bashforth’s method of 
experimenting with his chronograph and screens gives more trustworthy 
and more extensive information than the ballistic pendulum experiments 
of Hutton and Didion; and Hutton’s formula, as well as Didion’s, 
agrees with Mr. Bashforth’s Tables only for a limited range of velocities. 

Mr. Bashforth himself makes no attempt to condense his Tables into 
concise formula). Accordingly, as we shall presently see, he adopts a 
solution which is capable of being employed in conjunction with his 
Tables. He divides the trajectory into small arcs, and finds for every arc 
the time and the horizontal and vertical distances from one end of the 
arc to the other. The entire trajectory may then be plotted out, and the 
whole time and range may be discovered, as well as the final inclination 
of the direction of motion. The amount of labour, however, in calcu¬ 
lating all the quantities for a single component arc, even with the aid of 
copious tables, is so great that I was led to examine whether any thing 
could be done towards simplifying the solution and reducing the amount 
of calculation. It will appear in the sequel that rules of comparatively 
easy application can be employed, and that the tables necessary for their 
use are already existing or can be easily formed. 

Meanwhile, without entering into details, it will be convenient to give 
a brief account of the drift and scope of what is attempted in this paper. 

* Read May 11, 1876. See vol. xxy. p. 18. 
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The work is arranged under three heads, which are called the First, the 
Second, and the Third Methods, intended to signify three distinct and 
different solutions of the problem of finding the motion of shot. 

The First, which is the one adopted in Mr. Bashforth’s treatise, is a 
solution on the assumption that the retardation due to air-resistance is 
pv 3 , where fx is a constant. Now in the actual case fx is not constant; 
and therefore, in dealing with the equations of motion over any compo¬ 
nent part of the trajectory, a mean value of fx must be used. 

In contrast to the first method, the third adopts the mean value of 
another quantity, viz. the inclination of the direction of motion. The 
same thing was done in General Didion’s solution, in which it was neces¬ 
sary to employ the mean value of the cosine of the inclination. It will 
be found that in selecting this mean we really implicate more quantities 
than fx. But then there will be this advantage, that whereas there is no 
way of determining the mean of fx by the first method, and the greatest 
uncertainty prevails regarding it when comparatively large arcs are inte¬ 
grated over, according to the work now presented the nature of the mean 
is investigated. In fact the determination of the mean angle ' may 
be said to be the chief object and point in these investigations, because it 
will be seen that there is no difficulty in establishing any of the equations 
which will be used, and, excepting the labour, no difficulty in forming 
the Tables. 

What is described as the Second Method, although a distinct solution 
when the retardation is formulated by fxv n , is to be regarded as a mere 
stepping-stone to the one which follows it. It is in the second method 
that all the quantities are expressed in terms of the mean inclination and 
the magnitude of that angle determined. The quantity fx is taken con¬ 
stant ; and therefore, in the case when n = 3, if we use the same mean 
value of fx as in the first method, we ought to get a solution much the 
same. Beyond this, the second method possesses no further practical 
value in the business now in hand, being entirely superseded in that 
respect. Its chief value and importance consist in the determination of 
the mean angle, because it is shown that the same value of that mean may 
be used in the third method. 

This point being settled, it will be found that, according to the latter 
method, all the required quantities will depend on three integrals, two of 
which (the space- and time-integrals) have already been tabulated by Mr. 
Bashforth. The third, which may be called the velocity-integral, is now 
also tabulated for ogival-headed shot, and will be found further on. 

§ 2. It seems convenient at the outset to define the symbols which will 
be employed throughout the work :— 

v will denote the velocity of the shot at any point of the trajectory; 
u , the horizontal component of v ; 

the inclination of the direction of motion to the horizontal line in the 
plane of the trajectory (the deflection from the plane being neglected); 
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t, the time; 

w, the horizontal distance from some fixed point; 
y , the vertical distance. 

The integrations will be performed over a component arc of the trajec¬ 
tory, and the three last quantities are measured from the beginning of 
the arc. The values of t, oc, and y over the whole arc will be denoted by 
T, X, and Y. The values of u at the beginning and end of the arc will 
be denoted by^> and q; those of ty by a and /3. The acceleration due to 
gravity is denoted by g . 


Eirst Method. 

§3. The solution adopted by Mr. Bashforth is the famous one first 
given by John Bernoulli and published in 1721. It applies to any 
retardation formulated by yv n . All the characteristic quantities are 
expressed in terms of ty, which may be accomplished briefly thus :— 

The equations of motion in the horizontal direction and in the direc¬ 
tion of the normal to the trajectory in the ascending branch are 

~=- g v n co$<j>, .• (1) 


Hence 


city 

vJL = - gmH . 

. . . . (2) 

df 9 

. . . . (3) 

= iit(’‘+ 1 sec'‘ +1 4>. 

9 



Integrate this, recollecting that the initial values of u and ty are p and 
a, and get 

1 

u n 


1 if* 

«’* P" 9J<p 


n sec ,i+1 ty dty. 


P<P . , 

If the symbol P<p denotes 1 n sec n+l ty city, the equation will become 

Jo 

u n p g 


The horizontal velocity at the end of the arc is therefore given by 


i 

cf 



9 


. . (A) 
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At the vertex of the trajectory the velocity u 0 is given by 

JL = JL 4 . up • 

K p n 9 

« »;l </ ? / 


If we put 
we have 

From equation (2), 


_ n 

y v' 


(l — yY<p)n 


t= r a - 

J/» 9 


sec 


sec 2 ^ cty; 


by (5), 

Again, 


_ « 0 r 8ec2 ^# . 

= (i-yPo)** 


C?# 

dt 


\ a* 


Similarly, 


=J=~~ ~ S6c2 ^* 
9 J/3 ( 1 —yVp)» 

y_ f“ tan f sec 2 <jt d<j> 

9 J/3 (l — yP.fi)"* 


(4) 

(5) 


. (B) 


(C) 


(D) 


§ 4. The law of retardation found by Mr. Bashforth, in his own 
symbols, is expressed by 25i> 3 , where 


25 


d* K 
W (1000) 3 ’ 


d being the diameter of the cross section of the shot in inches, W its 
weight in pounds, and K a quantity which is a function of the velocity, 
and whose values are tabulated between certain values of the velocity for 
every ten feet. In integrating the equations of motion over any arc, the 
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quantity K, and therefore 2b, is taken constant, its mean value over the 
arc as near as it can be guessed at being used. 

Taking this law, the above results, written in the order in which they 
would be used by a calculator, are as follows 



. . (a) 

/ioooy_/ioooy_ & K p 

\ u 0 ) \ P ) W g 

• . (b) 

d- k/ w 0 y 

7 W g \1000/ ’. 

• • (c) 

m «■„ C sec 2 <l> d<i> 

<d(i-^) y ' ' 

• • (d) 

-y- l C a sec 

• • (e) 

u \ C a ^n <f) sec 2 (j> cty 

d Jp (i—r p <p) J 

• • (0 


The last three integrals Mr. Bashforth has given ample tables for, 
corresponding to different values of y and between certain ranges of 
angle. If there is no table for the exact value of y which results from 
equation (c), a method of interpolation must be employed. The integral 
P<p is in this case 3 tan <j> -f tan 3 <j>, and its values are tabulated as well as 
those of its logarithm. 

It will thus be seen that there are six distinct operations of some 
length, the first being the most serious, because there is some difficulty 
in getting K right. Supposing, however, that point to be settled—and I 
shall afterwards offer a few observations which will, I think, make the 
solution of (a) easier—the quantity y must be found. It will be seen that 
(b) and (c) are mere stepping-stones to the time- and distance-integrals. 

I shall now enter into an examination of the equations of motion, with 
the object of proving other formulae, which, when we have once discovered 
the velocity at the end of the arc, will give the time and distances with 
two operations. 

Second Method. 

§ 5. I remark, first, that although Bernoulli’s solution succeeds in 
expressing all the quantities in terms of integrals of 0, yet, owing to the 
difficulty and complexity of the integrals, it is practically valueless, except 
in two cases—first when n= 1, and next when w=3. I doubt, moreover, 
whether, in the case n — 3, it is the best solution when p is not constant, 
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and its mean value must therefore be taken over a small arc. The tables 
being such as they are, there can be no doubt that, if it were possible, the 
most convenient independent variable would be the velocity itself. On 
the other hand, all attempts at humouring the equations of motion so as 
to introduce the velocity as independent variable are of no avail. When 
the trajectory is very flat, it may be possible to get results which are not 
very objectionable; but no general solution is hereby attainable. The 
forms of the equations, however, suggest as a possibly good substitute 
the horizontal component of the velocity. Accordingly, taking this quan¬ 
tity as independent variable, I now proceed to find the distance-integrals 
by a method of approximation. 

Since 

dx 
dt 


and 




du 

dt 


= — fiv n cos (p. 


dx 1 n 

.\ — =--cos n 

du ji 


Similarly, 


.. X-lfcorf-V*. 

P J 2 u 

Y = i (* P sir 
w J t 


sin 6 cos 7 * 


du 


§ 6 . Our business now is to substitute for (j> its value in terms of u. 
To enable us to do this, put = and expand in powers of we 
have 

P a —P ? = ^_etc.; 


da 


_ _ JL = H (n sec “ +1 ad— sec " +1 a tan a^+ etc.Y ( 6 ) 

» jp n g \ 2 J 


Also 


cos 1 


i"- 1 0 =cos ”" 1 a + (n - 1 ) cos n ” 2 a sin a\p 


+ -—i. {(n— 2)cos w “ 3 aSin 2 a— C0S n l a}^ 2 . 
2 


>co 


4 -etc. 
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sin (f> cos' 4 ” 2 0=sin a cos' 4 ” 2 a + {(n~ 2) sin 2 a cos' 4 ” 3 a—cos' 4 ” 1 a^ 

-f {(w— 2)(n — 3) sin 2 a—(3w—5) cos 2 a} sin a cos' 4 ” 4 a . (8) 


I now propose to neglect the squares and higher powers of \p. The 
effect of this on the (c) integral will he that we have now to find X 
from 

cos' 4 " 1 a C p clu , (w —1) cos' 4 " 2 a sin a C p , du 

~i^ + -i=-i ^ 

Call the two integrals in this expression Q and It. Then 

^ n — 2^\cf~ 2 p n ~*) ’ 

and, taking account of equation (6), 

-o a cos"* 1 a C p / 1 1 1 V 7 


B=g £gg—g \ ( 

p n X v 

___g cos' 4 ^ 1 a f 1 
ji n \2n— 

Let this be put equal to 


Then, since 


we see that 


J 2 \M i> « V 

—2 \2 2n_2 j p 2n - 2 J n—2p u v 
equal to 

Q(«-P )/• 

_ JL\ 

t* n \t p"/ 


_1 _ 1 _ 

a ’* -2 


W 1 _ 

- 1 V- 1 1 f 

1 1 \ 

2n-2\q 2n ~ 2 

p2n-2j n — 2p U \ 

-Q I 
1 

^ ! 
8 

1 | 

_Lv 

'i_JL Vi_ 

■ A\ 

—2\ 

^«-2 p n ~ 2 ) \(f 

W 


Now let - — - =X, so that_pX^=^ —1'^ may be regarded as a frac- 
tion whose square may be neglected. We get 

« w +§^ w+ 7 ^-^:±^ w+ ... 

«(jA ) 2 /1 + (n - 2)p \+~ + 37 (p\f + . . \ 


n / •. n 2 i 3n 7^/ •> \3 * i 

|W+—g— 0*0 + 




l + («—2)pX+- 


— (i>X) 4 + ... 

W+ • • I’ 


1+VV 

,i l-K^-2) - x iiwe ° eglcct|yX),i 
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. f_i 3g+(3»—7)(_p— g) 

6 2+(»-2Xi>- 2 ) 

_j_(3»-7)p-(3«-10)? 

8 (w—2)p—(w—3)g ’ 

The value of X is now seen to become 

- {cos’*- 1 a+ (m-1) cos’*- 2 a sin «(«-/3)/}Q= icos’-^Q, 

fX fX 

where 

f—a—(a-(3)f 

_ a+/3. p — ( l U—fi 

2 ^(n—2}p-(«-3)2 # 6 * 

It is obvious that the value of Y may be obtained in exactly the same 
way, and that its magnitude is 

1 

-cos n “ 2 0 sin^Q. 

P 

§ 7. It is necessary now to enter into an examination of the magnitudes 
of the errors committed in neglecting squares of the quantity \p, and to 
discuss how far (that is to say, for what size of arc) we may with safety 
employ the values of X and Y just obtained. 

Before doing so, I shall briefly remark on the step that we took in 
neglecting (jp\) 2 in the value of/. On an examination of equation (9), it 
will be seen that the principal part of the error committed amounts to 

2n—5, 

—if 


The error, therefore, in the angle <j> is 


2n—5 (p — q) 2 
12 “ 




The consequent errors in X and Y, if we put D for the number of degrees 
in a— J3, are given by 


ax 

X 

aY 

Y 


= a similar expression. 




For the sake of simplicity, and to fix our ideas, let us take the case of 
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w=3, and suppose that we are integrating over an arc of 5°, in which one 
fourth of the horizontal velocity is lost. Then 

ax . - . 

_ — — tan (j> . ; 

*Y _ 

_=cot 2(p TF 1 ro - approximately. 

It will therefore be observed that the error committed is trifling. As 
<j> will be small over the greater part of any trajectory, this error will 
chiefly affect Y; but by good fortune the sign of the error is opposite 
to the error we shall presently find, and therefore helps to neutralize that 
error. 

Keturning to the integral for X, we are now to take account of the 
square of \p. Let the right-hand side of equation (6) be put equal to 

sec n+1 c vd/. 

9 

Then 

y=tf, tailed, 

A 

A 

Substituting this value of xp in the expansions of cos n_1 0 and cos n ~*(p sin <p, 
we get 

cos"” 1 a + (n — 1) cos n “ 2 a sin a xp f 

+ {2 n sin 2 a— -1} cos n “ 3 axp' 2 

A 

+ etc., 

and 

cos n ~ 2 a sin a + {(w—2) sin 2 a — cos 2 a} cos n ~ 3 axp 

-f {2(n—l)(w —2) sin 2 a — 4(w—1) cos 2 a} cos n-4 asin 

A 

+ etc. 

When the expressions just found are substituted in the (C) and (D) 
integrals, it is obvious that we have to determine these two integrals—- 
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Referring to the investigation of R in § 6 , we see that the work for R' 
will be the same, except that instead of a — /3 we should have to put an 
angle corresponding to it, such that 

i! = a —/3 - tan a(a-/3) 2 . 

The error in X due to the difference between R' and R is 

~l 

——— cos ”- 3 a sin 2 a(a — /3) 2 Q/. 

The integral S' may be reduced in the same way as the integral R. 
The most important part of it will be found to be 

^ 3 

The corresponding error in X is therefore 

^ (2n sin 2 a — 1) COS n ” 3 a Q. 

2fi o 

In the former of these two components of the error of X it will be 
sufficient for the purposes of an estimate to put/=|. In that case the 
sum of the two component errors (call it 5X) amounts to 

n-2 j-g n sin 2 «-1 _n +1 ^ a | co ^ a(a _ m 
= ^^(n-3)rinV-2| eOS **- 3 a(a-/3) 2 Q. 

It may be shown in a similar way that the error in Y is given by 
gY=^ (m ~ 1 )(w ~ 3 ) ; g naa ~ 5ffl + 11 J> sin a cos" _ 4 a(a—/3) 2 Q. 

A discussion of these expressions for any assigned value of n would 
determine for what magnitude of arc we might with safety employ the 
formulae for X and Y. I shall confine myself to the case when n = 3, and 
for the purposes of a ready estimate I shall take 

X= 1 cos 2 ^Q, 

v 1 a + /3 • a + /3 
x = - cos ~~~ sin — 

jj, 2 2 
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It thus appears that when w=3, 

SX_ 2 (a—/ 

X 14" cos (a ~f- /3) 3 

and 

SY__ 2 tang (a-ft) 2 
Y sin(a-f/3) 3 


Now suppose we were integrating over an arc of 5°, then 

o 

Light, approximately, be put equal to —§1-- =-l—. Our results would 


might, approximately, be put equal to —|A_ =jlp Our results would 
therefore be less than in error. Moreover the error in Y, which 
is really the more important of the two, is less than this, as I have 
pointed out at the beginning of this article. If, therefore, the formulae 
are otherwise serviceable, their inherent errors do not seem to be a great 
objection to their use. 

§8. The formulae for X and Y already found apply only to the 
ascending branch of the trajectory. A little consideration enables us to 
see that the same formulae apply to the descending branch, provided the 
mean angle </> is 


q + /3 _ P-~ C { ft—a 

2 (n— 2)%> — (n—3)q 6 * 

ft being now greater than a. 

§ 9. In the preceding articles we have neglected all consideration of 
the time-integral. It may, of course, be treated in the same way as the 
distance-integrals. I shall not, however, go into the general case, but 
shall merely state the result in the important case when w=3, 


nn_COS 2 ft C P clu 

where ft is equal to 

a ±fijL P~ ( l ( a ~ft) 

2 ‘ 9p — 3q 2 

in the ascending branch, and 

q+/3 __ p—q (ft—a) 

2 Qjp—Sq -2 

in the descending. 

If the arc integrated over is small, (ft is not very different from ft 
§ 10. I propose now to collect the results I have proved, stating them 
in the order and form in which they would be used in conjunction with 
Mr. Bashforth’s tables for K. 
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or 


or 


Summary of Results when n = 3, y— 


d 2 K 


W 100(f 


<rrJ-<rrf-$>-**■ 


2(j)=z a (ascending branch), 
P o 


a+/3 


(i-a 
P 8 


(descending branch), 


2^'=a+/3+ 42— (a—13) (ascending branch),. 

« + B — 0 > ~/^ (n — a) (descending branch); 
9p — 3q 


■v W,, , 0 -, 500,000/1000 1000\ 

X= ^( 1+ cos^)-^-(^— — 


JP 


1 • * 


w • 0 - 500,000/1000 1000\ 

Y=_ sm 2<b — ± —(---), . . . 

d K \ 2 P / 




!* _1200 

P 


2\ 


(a) 


(b) 


(o) 

(d) 

(e) 


These formulae might very easily be used if the calculator were furnished 
with tables of and f^^Y where N ranges between the magni- 


u; 


tudes of the velocity occurring in practice. 


Remarks on the Equation giving the Fall of Velocity in an Arc . 

§ 11. It will be observed that the two foregoing methods each open 
with the same equation (a). Now there is a serious difficulty in the use 
of that equation. Suppose, for example, we were to integrate over an arc 
of 1°: we should have to use the mean value of K between its values 
corresponding to the velocities at the beginning and end of the arc. But 
we do not know the latter of these velocities ; it is the very thing 
we have to find. The first steps in our work must therefore be to 
guess at it. The practised calculator can, from his experience, make a 
very good estimate. Having made his estimate, he determines K. He 
uses this value of K in equation (a); and if he gets the velocity he 
guessed at, he concludes that he guessed rightly and that he has got the 
velocity at the end of the arc. If the equation (a) does not agree with 
him, he makes another guess ; and so on, till he comes right. It seems 
to me, however, that this method of going to work, leaving out of account 
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the loss of time, is open to objection in the point of accuracy. For, first, 
there is no method of determining on what principle the mean value of 
K is to be found—what manner of mean it is. Again, let us suppose 
for an instant that the velocity at the end of the arc, guessed at, and the 
value of K are in agreement; that is to say, let the equation 

/lOOOV 8* /1000V 3 cl 2 K /T> , 

hold for the values of vp and K used by the calculator. It by no means 
follows that he has hit on the right vp and K. For if he is dealing with 

a part of the tables in which happens to be nearly equal to 


3W g sec 3 /3 (1000) 3 
F a Pj8 v 4 ’ 


it is obvious that there are ever so many pairs of values of vp and K 
which will stand the test of satisfying the above equation. Now an exa¬ 
mination of Mr. Bashforth’s tables for ogival-headed shot shows that the 
value of K diminishes as v increases from 1200 feet upwards, so that 
cZK 

• is negative for a considerable range of values of v which are common 

in practice. It is not at all unlikely, therefore, that the value for 
dK_ 

just stated may often be very nearly true; in which case the 

process of guessing becomes extremely dangerous. 

I shall in the next method give a plan for determining the velocity at 
the end of the arc, which seems to me simpler and more satisfactory than 
the one we have been now discussing. Meanwhile, in connexion with 
the present method, the following considerations are worthy of notice, 

§ 12. Let us consider the fundamental equation 

cU^d? K v 4 

d<p W g 1000| 3 

If ~ were uniform, it would be the fall of horizontal velocity in 

change of inclination equal to the unit of circular measure. As that 
angle is inconveniently large for discussions connected with trajectories 
of shot, let us take 1° for unit-angle, and let D be the number of degrees 
in the angle whose circular measure is <p. We then have 


W clu 

a? m j 


K 7T 
9 180 


x1000 x 



The quantity on the right-hand side is the same for all shot of the same 
kind; and I propose to find its values for intervals of the values of v. It 
will be convenient to represent those values by the ordinates of a curve 
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whose abscissae are the corresponding velocities. Besides the immediate 
purpose now in view, the same curve will be useful in exhibiting the 
manner of change of the retardation; since if y be the ordinate corre¬ 
sponding to velocity v, the corresponding retardation will be 

cl 2 180 gy 
W t r v 
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The primary use of these curves is that they enable us at a glance to 
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make a rough, estimate of the loss of velocity over an arc of one degree; 
and we have seen that such an estimate must first be made if we are to 
employ either of the foregoing methods. 

It is possible to make the preliminary guess more near the truth by 
considering the approximate character of the curves. [For example, the 
curve for ogival-headed shot is approximately a straight line, the tangent 
of whose inclination from 1700 down to 1250 is *5; and from 1250 down 
to 1000, *4; below 1000 it is unsafe to assign a value. "We can easily 
prove for fiat trajectories that the fall of velocity in one degree is either 

d 2 d 2 

W fclu\ (l-e W (du\ (l-e w \ 

-5 ) or -4 ) 

between the values above assigned. 


Third Method. 


§ 13. [Returning now to the fundamental equation 


clu Ev 



where E is the retardation due to the resistance of the air, since E is 
some function of v let us put it equal to f(y). Then, with our old nota¬ 
tion still in use, we get by integration 



clu 


-ft, 


clu 

u sec <pf(n sec (p ) 


ct—/5. 


Now instead of taking some mean value of the quantity K, as was 
done in the two previous methods, let us make the supposition that the 
quantity <p has its mean value—a supposition in this case by no means 
extravagant, since for the greater part of the trajectory sec <j> will vary 
very slowly. We then have, if D is the number of degrees in a —/3, 


180 cj rP 


1 


da 


q u sec <\>f(u sec (p ) 
Now put u=z Y co s </). The equation becomes 

sec<£ 




dV 


180? C 

- 3 q Sec <p W)' 


=1) sec 0. 
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With Mr. Bashf orth’s law substituted, this is 


d V d?- n - 

-ir (10 °0) _ Ky4=w Dsec ^ 

v or - 


• (a') 


1 q sec <J5 


Now the quantity 


180 <? 

7T 


iooofjl 


dY 

KY 4 


can be calculated for every 10 feet, by giving K its mean value over the 
10 feet, beginning with 1700 for ogival-headed shot and with 2150 for 
spherical. Supposing this were done, and the results for all the compo¬ 
nents of 10 feet added together, the number opposite any velocity v would, 
for example, in the ogival-head-shot tables, be the value of 


1800 , 
7 r 


/■»1700 


dY 

KY 4 * 


-ioOOf ( 

«/ V 

If we denote this quantity by V„, the equation (a') may be written 




sec <p ~ "Yp sec <p — ^ B Sec (p, • •••.. (a') 

which may be regarded as an equation for the determination of q. 

The distance-ordinates and time may be found in a similar manner. 
They are given by 


cP 

W 

d 2 

w 

cP 

w 


X = cos <j> 

f* 860 * 1000}W 
) - KY 2 ’ * * ' 

. » 

■ • (b’) 

« 

J q sec (p 



Y=sin $ 

r 8eo? lOOOfdV 
) ~ KY 2 ’ * * 

• • . 

• • (o') 

• 

J q sec <p 



T= 1 

c 

f sec * i000) 8 rfV 

KY 3 

1 q see <p 

. . . , 

• • (d’) 


If we denote by S v and T v the two integrals 


/•*1700 

J 

%/v 


umj\ s dY 

KY 2 


and 


✓*1700 

*Jv 


iooo] 3 c?v 

KY 3 ’ 


the equations (b'), (o'), and (d') become 


d 2 


yp X COS (j) (S ? sec <p Sp geo <p 


CO 
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cT 

w 

£ 

W 


Y — sin (j> (Bg. see <p —- Bp geo ), 

T—Tg seo<p “Tp S ec^* * 


* (C') 

* (d') 


Tlie quantities B v and T w have been tabulated by Mr. Bashforth for a 
considerable range of values of v, the upper limit being either 1700 
or 2150, according as the shot is ogival-headed or spherical. The 
tables for the ogival-headed shot have recently been revised and carried 
to one place further in decimals. The quantities B c and T e may 
therefore be fortunately regarded as completely determined; and the 
only question will be regarding the mean angle 0. Now it is a remarkable 
circumstance that the value obtained for $ in § 7, if q is not widely different 
from p, is very nearly the same for all values of n which are not far off from 
3. But for limited portions of the trajectory the retardation may be con¬ 
sidered as varying according to some simple power of the velocity, though 
that power is not the same from point to point, but still not far from 3, 
We may therefore take the value of <p found in § 7 as applicable to the 
method now in hand. The adoption of this value of the mean angle, 
since Y=Xtan <j>, is really equivalent to supposing the shot to move 
parallel to the chord; and the above proof shows what the limits of inte¬ 
gration must be in order that the supposition may be made to approximate 
to the actual case. The most sensitive quantity in this method, especially 
near the vertex of the trajectory, is Y; and in finding it over an arc cor¬ 
responding to a change of inclination as large as 5°, it is necessary to use 
the correct value obtained in § 7. As regards the other three quantities 
given by the integrals (a'), (b'), and (d'), it will not matter much if w r e 


take 


a + /3 
2 


, at least for flat trajectories. 


If, however, the trajectory is not 


flat, and extreme accuracy is needful, it will be necessary to determine by 
(a') the quantity q twice over—first an approximate value of it, in order 
to get the mean angle next its value with the mean angle used in the 
equation. 

The following Table gives the value of V„ for ogival-headed shot as low 
as 900 feet per second. Below that value of the velocity Mr. Bashforth 
does not give tables for the resistance, and the magnitudes of the resist¬ 
ance for the lower velocities of ogival-headed shot have yet to be found. 
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Table of Values of Y v = 


180 9 

TV 



dV 

KV 4 " 



0 . 

1 . 

2 . 

3 . 

4 . 

5 . 

6 . 

7 . 

8 . 

9 . 

90 .. . 

91 .. . 

92 .. . 

93 .. . 

94 .. . 

95 .. . 

96 .. . 

97 .. . 

98 .. . 

99 .. . 
100 ... 
101 ... 
102 ... 

103 .. . 

104 .. . 

105 .. . 

106 .. . 

107 .. . 

108 .. . 

109 .. . 

no ... 

111 .. . 
112 ... 

113 .. . 

114 .. . 

115 .. . 

116 .. . 

117 .. . 

118 . J 

119 .. . 

120 .. . 
121 ... 
122 ... 

123 .. . 

124 .. . 

125 .. . 

126 .. . 

127 .. . 

128 .. . 

129 .. . 

130 .. . 

131 .. . 

132 .. . 

133 .. . 

134 .. . 

135 .. . 

136 .. . 

137 .. . 

138 .. . 

139 .. . 

140 .. . 

141 .. . 

142 .. . 

143 .. . 

144 .. . 

145 .. . 

8*3380 

7*9125 

7-5006 

71213 

6-7596 

6-4146 

6-0880 

5-7797 

5-4898 

5-2169 

4-9623 

4-7252 

4-5050 

4-3016 

4-1148 

3-9439 

3-7871 

3*6422 

3-5069 

3-3792 

3-2576 

3-1413 

3-0296 

2-9223 

2-8191 

2-7198 

2-6244 

2-5319 

2-4429 

2-3570 

2-2739 

2-1936 

2-1159 

2-0406 

1-9678 

1-8972 

1-8287 

1-7625 

1-6981 

1-6356 

1-5749 

1-5157 

1*4584 

1-4026 

1*3484 

1 * 295 £ 

1 - 244 C 

1 T 93 £ 

1 * 144 £ 

1 - 097 C 

l * 050 c 

1-0047 

•9603 

* 916 £ 

•8737 

*8311 

8-2955 
7-8714 
7-4627 
7-0851 
6-7251 
6-3819 
6*0571 
5-7507 
5-4625 
5-1914 
4-9386 
4*7032 
4-4847 
4*2829 
4-0977 
3-9282 
3*7726 
3-6287 
3-4941 
3-3670 
3-2460 
31301 
3-0189 
2-9120 
2*8092 
2-7102 
2-6152 
2-5230 
2-4343 
2-3486 
2-2659 
2-1858 
2-1084 
20333 
1-9607 
l 1-8903 
1-8221 

> 1-7560 
1-6919 

■ 1*6295 
i 1-5689 
' 1*5099 
: 1 * 452 £ 
1 1-3972 
t 1-3431 

> 1-2904 
) 1 - 239 C 
5 1*1882 
5 1 - 140 C 
) 1 - 092 E 
5 1-0457 
r 1*0002 
. -9557 

> -9122 
r - 869 £ 
) - 827 £ 

8-2529 
7-8302 
7-4248 
7-0489 
6-6906 
6*3492 
6*0263 
5-7217 
5-4352 
5*1660 
4*9149 
4*6812 
4-4643 
4*2642 
4-0806 
3-9125 
3-7581 
3-6151 
3*4814 
3-3549 
3-2344 
3-1190 
3-0081 
2*9017 
2-7992 
2-7007 
2-6059 
2-5141 
2-4257 
2-3403 
2-2578 
2-1781 
2-1008 
2-0260 
1-9537 
1-8834 
1-8154 
• 1-7496 
i 1-6856 
. 1-6234 
i 1-5630 
l 1-5042 
; 1*4472 
1 1*3915 
. 1-3375 
: 1-2852 
) 1 - 234 C 
) 1 - 184 C 
) 1-1352 
\ 1-0876 
r 1-0412 
J *9955 
r -9514 
5 -9072 
5 -8655 
5 -8237 

8-2104 
7-7890 
7*3868 
7*0127 
6-6561 
6-3166 
5-9954 
5-6927 
5-4079 
5-1405 
4-8912 
4-6591 
4-4440 
4*2456 
4-0635 
3*8968 
3*7436 
3-6016 
3*4686 
3-3427 
3-2228 
3*1078 
2-9974 
2-8913 
2-7893 
2-6916 
2*5967 
2-5052 
2-4171 
2-3320 
2-2498 
21703 
2-0933 
2*0188 
1-9466 
: 1-8766 
: 1-8088 

> 1-7432 
■ 1-6794 
: 1-6174 
* 1-5571 
! 1-4984 
! 1-4417 
S 1-3863 
S 1*3325 
! 1-2801 
) 1-2289 

> 1*1791 
! 1-1305 
i 1-0829 
5 1-0366 
1 -9915 
t - 947 C 
) -9036 
5 -8619 
r -8196 

8-1678 
7-7479 
7*3489 
6*9766 
6*6216 
6-2839 
5‘9646 
5*6637 
5*3806 
5-1151 
4*8675 
4*6371 
4-4236 
4-2269 
4-0464 
3-8811 
3-7292 
3-5881 
3-4558 
3-3306 
3*2111 
3-0966 
2-9867 
2-8810 
2-7794 
2-6816 
2-5874 
2-4963 
2-4085 
2-3237 
2*2418 
2-1625 
2-0858 
2-0115 
1-9396 
1-8697 
1-8022 
1-7367 
1-6731 
1-6113 
1*5512 
: 1*4927 
1-4361 
i 1-3809 

> 1-3272 
. 1-2749 

> 1*2239 
1-1742 

> 1*1257 
) 1*0782 
; 1*0321 
; -9869 
) -9427 
5 *8994 
5 * 857 C 
» -8155 

81253 
7*7066 
7*3110 
6-9404 
6*5870 
6*2513 
5*9338 
5*6347 
5*3533 
5*0896 
4-8437 
4*6151 
4*4033 
4*2082 
4*0293 
3*8655 
3*7147 
3-5746 
3-4431 
3-3184 
3-1995 
3-0855 
2-9760 
2-8707 
2*7695 
2-6721 
2-5782 
2-4874 
2-3999 
2*3154 
2-2338 
2*1548 
2*0783 
2*0042 
1*9325 
1*8629 
1*7956 
1-7303 
1*6669 
1*6052 
! 1-5453 
1*4870 
1*4305 
i 1*3755 
! 1*3220 
i 1*2698 
| 1-2189 
! 1*1692 
r 1*1209 
; 1-0736 
1*0275 

> -9824 
' *9382 
: -8951 
) *8522 

► -sin 

8-0827 
7-6654 
7-2730 
6-9043 
6*5525 
6-2186 
5*9030 
5-6057 
5-3260 
5-0641 
4-8200 
4*5931 
4*3830 
4*1895 
4-0122 
3-8498 
3-7002 
3-5610 
3*4303 
3-3062 
3-1879 
3*0743 
2-9652 
2-8604 
2-7595 
2-6625 
2-5689 
2-4785 
2-3913 
2-3071 
2-2257 
2-1470 
2-0707 
1-9969 
1-9254 
1-8561 
1-7889 
1-7239 
1-6606 
1-5991 
1-5394 
1*4813 
1*4249 
1-3701 
1*3167 
1*2646 
1*2139 
T 1644 
1-1161 
1-0689 
1-0229 
•9779 
•9339 
. -8908 
\ -8488 
t * 807 £ 

8-0402 
7-6242 
7*2351 
6 * 868 ] 
6-5180 
6*1860 
5*8722 
5-5767 
5-2987 
5*0387 
4-7963 
4-5711 
4-3626 
4-1708 
3*9951 
3*8341 
3-6857 
3*5475 
3*4175 
3-2941 
3-1762 
30631 
2*9545 
2*8501 
2*7496 
2-6530 
2-5597 
2-4696 
2*3827 
2*2988 
2*2177 
2*1392 
2*0632 
1-9896 
1-9184 
1*8492 
1*7823 
1-7174 
1*6544 
1*5931 
1-5334 
i 1-4756 

* 1-4193 
1*3646 
1-3114 

> l - 259 £ 

• 1-2089 
: 1 * 159 £ 
. 1-1113 

> 1-0643 
) 1*0184 
) *9733 

> -9296 
J -8866 

> -8444 
! *8035 

7*9976 
7-5830 
7*1972 
6-8319 
6-4835 
6-1533 
5*8413 
5-5477 
5-2715 
5*0132 
4*7726 
4-5490 
4*3423 
4-1522 
3-9780 
3*8184 
3*6712 
3-5340 
3-4047 
3-2819 
3-1646 
3*0519 
2-9438 
2-8397 
2-7397 
2-6435 
2-5504 
2-4607 
2-3741 
2-2905 
2-2097 
2-1314 
2-0557 
1*9824 
1-9113 
1-8424 
1*7757 
1-7110 
1*6481 
1*5870 
: 1*5275 

> 1*4698 
; 1-4138 

> 1*3592 
: 1-3061 
* 1*2543 
i 1-2038 

> 1*1546 
l 1-1066 
; 1*0596 
: 1-0138 
; - 969 C 
; -9252 

> *8823 
L *8402 
5 -7991 

7*9551 
7-5418 
7-1592 
6-7958 
6-4490 
6*1207 
5-8105 
5-5187 
5-2442 
4-9878 
4-7489 
4-5270 
4*3219 
4*1335 
3*9609 
3-8027 
3-6567 
3-5204 
3-3920 
3-2698 
3-1529 
3*0408 
2*9330 
2-8294 
2*7297 
2-6339 
2-5412 
2-4518 
2-3655 
2*2822 
2*2016 
2-1237 
2*0481 
1-9751 
1-9043 
1-8355 
1-7690 
1*7045 
1*6419 
1*5810 
1*5216 
1-4641 
1*4082 
1-3538 
1-3008 
1*2492 
1-1988 
1-1497 
1-1018 
1*0549 
1*0093 
•9646 
•9209 
•8780 
•8361 
. -7950 
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Table ( continued ). 



0 . 

1 . 

2 . 

3 . 

4 . 

5 . 

6 . 

7 . 

8 . 

9 . 

146 ... 

•7910 

•7870 

•7830 

•7790 

•7750 

•7710 

•7669 

•7629 

•7589 

•7549 

147 ... 

•7509 

•7470 

•7430 

•7391 

•7352 

•7312 

•7273 

•7234 

•7194 

•7155 

148 ... 

•7115 

•7077 

•7038 

•6999 

•6961 

•6922 

•6884 

•6845 

•6806 

•6768 

149 ... 

•6729 

•6691 

•6654 

•6616 

•6578 

•6540 

•6502 

•6465 

•6427 

•6389 

150 ... 

•6351 

•6314 

•6277 

•6239 

•6202 

‘6165 

•6128 

•6091 

•6053 

•6016 

151 ... 

•5979 

•5943 

•5906 

•5870 

*5833 

•5797 

•5760 

•5724 

•5687 

•5651 

152 ... 

•5614 

•5578 

•5542 

•5506 

•5470 

•5435 

•5399 

•5363 

•5327 

•5291 

153 ... 

•5255 

•5220 

•5185 

•5149 

•5114 

■5079 

•5044 

•5009 

•4973 

•4938 

154 ... 

•4903 

*4868 

•4834 

•4799 

•4765 

•4730 

•4695 

•4661 

•4626 

•4592 

155 ... 

•4557 

•4523 

•4489 

•4455 

•4421 

•4387 

•4352 

•4318 

•4284 

•4250 

156 ... 

•4216 

•4183 

•4149 

•4116 

•4082 

•4049 

•4015 

•3982 

•3948 

•3915 

157 ... 

•3881 

•3848 

•4815 

•3782 

•3749 

•3716 

•3684 

•3651 

•3618 

•3585 

158 ... 

•3552 

•3520 

•3487 

•3455 

•3422 

•3390 

•3357 

! -3325 

•3292 

•3260 

159 ... 

•3227 

•3195 

•3163 

•3132 

•3100 

•3068 

•3036 

•3004 

•2973 

•2941 

160 ... 

•2909 

•2877 

•2846 

•2814 

•2783 

•2751 

•2719 

1 -2688 

•2656 

•2625 

161 ... 

•2593 

•2562 

•2531 

•2500 

•2469 

! *2439 

•2408 

•2377 

•2346 

•2315 

162 ... 

•2284 

•2254 

•2223 

•2193 

•2162 

! *2132 

•2102 

•2071 

•2041 

•2010 

163 ... 

•1980 

• 1950 ; 

•1920 

•1890 

•1860 

! T831 

•1801 

•1771 

•1741 

•1711 

164 ... 

•1681 

•1652 

•1622 

•1593 

•1563 

•1534 

•1505 

•1475 

•1446 

•1416 

165 ... 

•1387 

•1358 

•1329 

•1300 

•1271 

•1243 

•1214 

•1185 

•1156 

•1127 

166 ... 

•1098 

•1070 

•1042 

•1013 

•0985 

•0957 

•0929 

•0900 

•0872 

•0844 

167 ... 

•0816 

•0788 

•0760 

•0732 

•0705 

♦0677 

•0649 

•0621 

•0594 

•0566 

168 ... 

•0538 

! - 0511 ; 

•0484 

•0456 

•0429 

•0402 

• 0375 ; 

•0347 

•0320 

•0293 

169 ... 

•0266 

•0239 

•0213 

•0186 

•0159 

•0133 

•0106 

•0080 

•0053 

•0027 


An example of the use of these Tables was given in the abstract which 
was printed in the e Proceedings ’ of the Society, vol. xxv. p. 18. 


PosrscKtrTe 

Professor «T. Couch Adams, to whom this paper was shown before 
publication, has obtained a solution of the equations employed in the 
second method, which is of a more complete and satisfactory character 
than the one given above. The results he has arrived at are contained 
in the following note, which he has kindly allowed me to subjoin :~— 

“ Employing the notation of the paper, and supposing the resistance to 
vary as the wth power of the velocity, the horizontal velocity q is given 
by the equation 

I (sec^J +I («~/3){1+AO + !) [OH-2) (sec^ J-(n+l)j(a-/3) 2 

where a —/3 is expressed in the circular measure. 

The inclination </> of the chord AB is given by 

where a— $ is supposed to be expressed as before in the circular measure. 

If a—jo be expressed in minutes, the last term must be multiplied by 
sin Ih 
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M The value of the mean angle <f> to be employed in finding the time- 
integral is 



where the last term is the same as that in the above value of </>, but the 
second term is only one half of its amount in the former case. 

“ It will be seen that the above expressions for $ and <j> are independent 
of the value of n* 

“ Also, if 

Q=Wt(«- 1) [(»-2) (sec a + /5 ) - («~3)](«-/3) 2 , 

where a—ft is expressed as before in the circular measure, the values of 
the coordinates X, Y, and of the time T are given by 

T= -..—.l. Q(cos 0 ' )’— 1 ( ~— iA 

VI P'J 

14 It may be remarked that if a and ft, as well as p and q, be interchanged, 
the values of <£, <j>, and Q will remain unaltered, and the values of X, Y ? 
and T will merely change their signs, as it is evident should be the case. 

44 The above values of 7p, <p\ and Q are true to the third order of small 
quantities inclusive, and the values of X, Y, and T are true to the fourth 

order, considering and a—ft to be small quantities of the first 

order/’ 


von. xxvi. 
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